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Wave functions and form factors of vector mesons are investigated in the holographic dual model 
of QCD with oscillator-like infrared cutoff. We introduce wave functions conjugate to solutions 
of the 5D equation of motion and develop a formalism based on these wave functions, which are 
very similar to those of a quantum-mechanical oscillator. For the lowest bound state (p-meson), we 
show that, in this model, the basic elastic form factor exhibits the perfect vector meson dominance, 
i.e., it is given by the p-pole contribution alone. The electric radius of the p-meson is calculated, 
{r 2 p )c = 0.655 fm 2 , which is larger than in case of the hard- wall cutoff. The squared radii of higher 
excited states are found to increase logarithmically rather than linearly with the radial excitation 
number. We calculate the coupling constant f p and find that the experimental value is closer to 
that calculated in the hard-wall model. 
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I. INTRODUCTION 

Holographic duals of quantum chromodynamics 
(QCD) which are based on the gauge/gravity correspon- 
dence |l| have been applied recently to hadronic physics 
(see, e.g., I I II, 0, i, I M EL M M M), 
and demonstrated their ability to incorporate such es- 
sential properties of QCD as confinement and chiral sym- 
metry breaking, and have been successful in many cases 
in determination of static hadronic properties, i.e., reso- 
nance masses, decay constants, chiral coefficients, etc. 
In Refs. @, H|, the dynamic properties (form factors) 
have been studied within the holographic approach of 
Ref. and the connection between AdS/QCD approach 
of Refs. @, [|[ and the usual light-cone formalism for 
hadronic form factors was proposed in [TTI | and discussed 
in [l5| . The calculation of form factors of scalar and vec- 
tor hadrons within the approach of Ref. 0] was performed 
in Refs. [H, [I3|, and applied to study the universality of 
the p-meson couplings to other hadrons. In our recent 
paper fl8| , we studied form factors and wave functions of 
vector mesons within the framework of the holographic 
QCD model described in Refs. @, 0, H] (which will be 
referred to as hard- wall model). 

In the hard- wall model , the confinement is modeled by 
hard- wall cutting off the AdS space along the extra fifth 
dimension at some finite value z = Zq. The solutions of 
the relevant eigenvalue equation are given by the Bessel 
functions, and masses of bound states are given by the 
roots M n = jo,n/zo of Jq{Mzo). As a result, the masses 
of higher excitations behave like M 2 ~ n 2 . It was argued 
[l2[ [l9j that, instead, one should expect M 2 ~ n behav- 
ior. This connection can be derived from semiclassical 
arguments [l!|[2(|. An explicit AdS/QCD model which 
gives such a linear behavior was proposed in Ref. 
The hard-wall boundary conditions in this model are sub- 
stituted by an oscillator-type potential providing a soft 
IR cut-off in the action integral (for this reason, it will 



be referred to as "soft- wall model"). 

The aim of the present paper is to study form factors 
and wave functions of vector mesons within the frame- 
work of the soft-wall model formulated in Ref. . and 
compare the results with those we obtained in Ref. [18| in- 
vestigating the hard-wall model . To this end, we extend 
the approach developed in Ref. [l8| . We start with recall- 
ing, in Section II, the basics of the soft-wall model and 
some results obtained in Ref. [l~2fl . in particular, the form 
of the relevant action, the eigenvalue equation for bound 
states and its solution. In Section III, we derive a useful 
integral representation for the bulk-to-boundary propa- 
gator V(p, z) that allows to write V(p, z) as an explicit ex- 
pansion over bound state poles with the z-dependence of 
each pole contribution given by ll tp wave functions" that 
are eigenfunctions of the 5D equation of motion. Then we 
show that the same representation can be obtained from 
the general formalism of Green's functions. However, 
as we already emphasized in Ref. the ip n (z) wave 
functions are not direct analogues of the usual quantum- 
mechanical wave functions. In particular, a meson cou- 
pling constant f n is obtained from the derivative olij) n (z) 
at z = rather than from its value at this point. To this 
end, we introduce wave functions" which look more 
like wave functions of oscillator bound states in quantum 
mechanics. Their values at z = give the bound state 
couplings g^fn/Mn, they exponentially decrease with z 2 , 
and thus they have properties necessary for the light-cone 
interpretation of AdS/QCD results proposed in Ref. [TT[ | . 
In Section IV, we study the three-point function (J J J) 
and obtain expression for transition form factors that in- 
volves ip wave functions and the nonnormalizable mode 
factor Sf(Q,z). The latter is written as a sum over all 
bound states in the channel of electromagnetic current, 
which gives an analogue of generalized vector meson dom- 
inance (VMD) representation for hadronic form factors. 
In Section V, it is shown that it is possible to rewrite 
form factors in terms of 6 functions. Then we formu- 
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late predictions for p-meson form factors, and analyze 
these predictions in the regions of small and large Q 2 . 
In particular, our formalism allows to calculate p- meson 
electric radius, and the radii of higher excited states. It is 
also shown that, for the basic p-meson form factor T(Q 2 ) 
given by the overlap of the <f> wave functions, the soft-wall 
model predicts exact VMD pattern, when just one lowest 
bound state in the Q 2 -channel contributes. For another 
p-meson form factor F(Q 2 ), which is given by the over- 
lap of the ip wave functions, a two-resonance dominance 
is established, with only two lowest bound states in the 
Q 2 -channel contributing. In Section VI, we compare our 
results obtained in the soft-wall model with those derived 
in the hard- wall model studies performed in Ref . [HI . At 
the end, we summarize the paper. 



II. PRELIMINARIES 

We consider the gravity back grou nd with a smooth 
cutoff that was proposed in Ref. [12| instead of a hard- 
wall infrared (IR) cutoff. In this case, the only back- 
ground fields are dilaton x( z ) — z 2 n 2 and metric gMN- 
The metric can be written as 



MjN 



gMNdx dx 



(rj fiv dx )1 dx 1 ' — dz 2 



(1) 



where r]^ — Diag(l, — 1, — 1, — 1) and p, v = (0,1,2,3), 
M, N = (0, 1, 2, 3, z). To determine the spectrum of vec- 
tor mesons, one needs the quadratic part of the action 



S 



1 



AdS 



4 dz _ 
ax — e 

z 



Tr (F MN F 



MN\ 



(2) 



where F M n = d M V N - d N V M - i[V M , V N ], V M = t a V^, 
(t a = <j a /2, with a a being Pauli matrices). In the axial- 
like gauge V z = 0, the vector field VHx,z = 0) corre- 
sponds to the source for the vector current J®(x). To 
obtain the equations of motion for the transverse compo- 
nent of the field, it is convenient to work with the Fourier 
transform V"(p, z) of V®(x, z), for which one has 



d z 



e~ z d z K(p,z) 



p 2 l e - z2 V«(p,z) 







(3) 



(Here, and in the rest of the paper, we find it convenient 
to follow the convention of Ref. [l2j , in which the oscilla- 
tor scale k is treated as 1, i.e., we write below z 2 instead 
of k 2 z 2 , e~ z instead of e~ z K , etc. Using dimensional 
analysis, the reader can easily restore the hidden factors 
of k in our expressions. In some cases, when n is not 
accompanied by z, we restore k explicitly.) The eigen- 
value equation for wave functions ipn(z) of the normal- 
izable modes can be obtained from Eq. by requiring 



P 



M 2 , which gives 



l 



Z tpn 



~4>n = 



(4) 



As noted in Ref. [T2|, the substitution 

1> n (z) = e z2 / 2 ^z~y n (z) (5) 
gives a Schrodinger equation 

•! \ , . ., . (6) 



K + [ z + ^2 ) *« = M n *n , 



which happens to be exactly solvable. The resulting spec- 
trum is M 2 = 4(n + 1) (with n = 0,1,...), and the 
solutions ip n (z) of the original equation ^ are given by 



ipn(z) = z' 



n+1 



Li(z 2 ) , 



(7) 



where L x n (z 2 ) are Laguerre polynomials. The functions 
ip n (z) are normalized according to 



— e z2 ip m (z) ijj n (z) =6„ 



(8) 



Correspondingly, the ^ n (z) functions of the Schrodinger 
equation |6j) arc normalized by 



dz^ m {z) *n(z) = Smn , 



(9) 



i.e., just like wave functions of bound states in quan- 
tum mechanics. Note, however, that the functions ^ n (z) 
behave like z 3 ^ 2 for small z, while quantum-mechanical 
wave functions of bound states with zero angular momen- 
tum have finite non-zero values at the origin. 



III. BULK-TO-BOUNDARY PROPAGATOR 

It is convenient to represent V£(p, z) as the product 

of the 4-dimensional boundary field (p) and the bulk- 
to-boundary propagator V(p,z) which obeys the basic 
equation 



d z 



z d z V 



2 1 - 
p — e 

z 



V = 



(10) 



that follows from Eq. j3| and satisfies the boundary con- 
dition 



V( P> z = 0) = l 



(11) 



Its general solution is given by the confluent hypergeo- 
mctric functions of the first and second kind 



V(p,z) = A ^{a^.z^+B f7(a,0,z 2 ) , 



(12) 



where a = — p 2 /4k 2 , A and B are constants. Since the 
function \Fi(a, 0, z 2 ) is singular for z = 0, we take A = 0. 
Then, for a > 0, the bulk-to-boundary propagator V{p, z) 
can be written as 



V(p, z) — a I dxx" 
'o 



exp 



(13) 
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It is easy to check that this expression satisfies Eqs. ([TO]) 
and Integrating by parts produces the representa- 

tion 



V(p,z) = z 2 



dx 



■ x exp 



1-x 



> (14) 



from which it follows that if p 2 — (or a = 0), then 

V(0,z) = l (15) 



for all z. The integrand of Eq. (114]) contains the gener- 
ating function 



1 



{l-xf 



exp 



l-at 



5> r \(z 2 )x" (16) 



n=0 

1 



for the Laguerre polynomials L n {z ), which gives the rep- 
resentation 



n=0 



a + n + 1 



(17) 



that can be analytically continued into the timelike a < 
region. One can see that V(p, z) has poles there at ex- 
pected locations p 2 = 4(n + 1)k 2 . 

The same representation for V(p, z) can be obtained 
from the Green's function 



G(p;z,z') = J2 



1pn{z)^n{z') 

p 2 - Ml 



corresponding to Eq. (|10p . namely, 
1 



V(p,z') 



-e z d z G(p; z, z') 



(18) 



(19) 



E 



^ P 2 -M 2 ' 



P 2 - Ml 



n=0 



which coincides with Eq. (| 1 7[) . 

The two-point density function can also be obtained 
from the Green's function: 



5V) 



1 



E 



i 



e- z "a z G(p;z,z') 



f 2 

J n 



^ n P 2 -M 2 n ' 



(20) 



n=0 



where the coupling constants /„ = K 2 i/8(n + l)/<?5 ob- 
tained in [12j are determined by 



fn= e Z d z tp n (z) 

The propagator V(p, z) can be represented now as 



V{p,z) = 35 E 



fnlpn(z) 

M 2 -p 2 



(21) 



(22) 



where ip n (z) are the original wave functions ([7]) corre- 
sponding to the solutions of the eigenvalue equation 

Given the structure of Eq. (|2"Tj) , it is natural to intro- 
duce the conjugate wave functions 



(f>n{z) 



M„z 



e z d z ip n (z) 



^-^(z 2 )-^ 2 ^ 2 )] , (23) 



whose nonzero values at the origin f n 95/M n are propor- 
tional to the coupling constant /„ (in this particular case, 
fn9b/M n = s/2k). The inverse relation between the ip 
and d> wave functions 



1 2 

4>n{z) = ~ — ze z d z <p n {z) 
IVl n 



(24) 



can be obtained from Eq. (U). The 0- functions are nor- 
malized by 



dzze Z (f> m (z) 4>n{z) = 5 V 



(25) 



In particular, for the lowest states, we have 

o (z) = V2e" z2 , ^(z) = V^e^ 2 (1 - z 2 ) . (26) 

Just like zero angular momentum oscillator wave func- 
tions in quantum mechanics, these functions have finite 
values at z = 0. They also have a Gaussian fall-off er z 
for large z. To make a more close analogy with the oscil- 
lator wave functions, it makes sense to absorb the weight 
e z in Eq. (|25[) into the wave functions, i.e., to introduce 
"<&" wave functions 

$„(«) = e z2 / 2 ^ n (z) = -}—e- z "l 2 d z i, n {z) , (27) 
M n z 

which are nonzero at z = 0, decrease like e^ 2 ^ I 2 for large 
z, and are normalized according to 



f 

Jo 



dzz$ m (z) $ n (z) = 8 n 



(28) 



The presence of the z weight in this condition (which 
cannot be absorbed into wave functions without spoiling 
their behavior at z = 0) suggests that pursuing the anal- 
ogy with quantum mechanics one should treat z as the 
radial variable of a 2-dimensional quantum mechanical 
system. 



IV. 3-POINT FUNCTION 



The variation of the trilinear (in V) term of the action 
S { a1s = ~ e -§J § ^ (W) V» h V» (29) 



4 



calculated on the solutions of the basic equation (fT0|) 
gives the following result for the 3-point correlator: 

{JaiPi)^{-P2)^{q)) = e abc ^) A ^5^{ Pl -p 2 + q ) 
xT a ^( Pl ,p 2 ,q)W(p uP2 ,q) , (30) 
with the dynamical part given by 

W{pi,p2,q)= / — e- z2 V(p 1 ,z)V(p 2 ,z)V(q,z) , 

(31) 

and the kinematical factor having the structure of a non- 
abelian three-field vertex: 

T a ^( Pl , P2 , q) = V a "(Q ~ Pif ~ V^(P2 + q) a 

+ 7] aP (Pi+P2Y • (32) 

Incorporating the representation Eq. (j22f for the bulk- 
to-boundary propagators gives the expression 

rrl 2 2 r^i2\ fnfkF n k(Q ) / 00 \ 

T( Pl ,p 2> Q)= E i(p? _ M »)(pS- Mi) 

for T(p\,p 2 ,Q 2 ) = W(pi,p2,q)/gl as a sum over the 
poles of the bound states in the initial and final states. In 
the z-integral of Eq. (f3"Tj) . the contribution of each bound 
state is accompanied by its wave function ip n {z), while 
the g-channel is represented by J(Q, z) = V(iQ, z). This 
gives the Q 2 -dependent coefficients 

F nk (Q 2 )= —e- z2 J(Q,z)Mz)Mz), (34) 
Jo z 

which have the meaning of transition form factors. Note 
that since J~(0,z) = 1, the orthonormality relation © 
assures that F nn (Q 2 = 0) = 1 for diagonal transitions 
and F nk (Q 2 = 0) = if n ^ k. 

The factor CT{Q, z) can be written as a sum of 
monopole contributions from the infinite tower of vector 
mesons: 



J(Q,z)=9sJ2 



frnlprn(z) 

Q 2 + Ml 



(35) 



This decomposition, discussed in Refs. [16|, LL8|, directly 
follows from Eq. (f2"2"| . As a result, the form factors 
F nk (Q 2 ) can be written in the form of a generalized VMD 
representation: 



Fnk (Q 2 ) = E I 



i 1 + Q 2 l M m ' 



(36) 



where the coefficients F m ^ nk are given by the overlap in- 
tegrals 

F m ,nk = / — e~ z tp m {z) , ipn{z) , ipk{z) ■ (37) 

JH m .In Z 



V. FORM FACTORS 

In terms of the 'J wave functions of the Schrodinger 
equation ([6]), the form factors are given by 



F n k(Q 2 ) 



dzJ(Q,z)V n (z)V k (z) 



(38) 



which looks like an expression for form factors in quan- 
tum mechanics. However, as we discussed above, the '3/ 
wave functions are not direct analogues of quantum me- 
chanical wave functions. For such an analogy, the $ wave 
functions l[27p are much more suitable objects. So, let us 
introduce form factors involving $ wave functions 



Fnk{Q 2 ) 



dzzJ{Q,z) $„(z)$ fc (z) 



(39) 



Again, since J7(Q = 0, z) = 1 for all z, the normalization 
condition ([28]) for the & n (z) wave functions guarantees 
that the diagonal form factors ^F n n{Q 2 ) are normalized 
to 1 for Q 2 — 0, while the non-diagonal ones vanish when 
Q 2 = 0. To establish connection with F nk (Q 2 ) form fac- 
tors, we use Eq. (|27[) to substitute $ functions by deriva- 
tives of ip wave functions, which gives 

/>00 7 

M n M k T nk {Q 2 ) = \ — e-* 2 J(Q, z) tf n {z) ^' k (z) . 



(40) 

Integrating ip' k by parts, taking into account that 
(0) = and incorporating the eigenvalue equation (|3|) 
for ip n gives 



M n M k T nk {Q 2 ) = M 2 F nk (Q 2 ) 



(41) 



-e~ z iP' n (z)Mz)d z J(Q,z) 



Similarly, integrating tp' n by parts we obtain 
M n M k T nk {Q 2 ) = M 2 F nk (Q 2 ) (42) 
dz 



— e- z Mz)*P'k(z)d z J{Q,z) . 



Adding these two expressions, integrating (ip n ipk)' by 
parts and using the basic equation (JTOJ) for J(Q, z) gives 



F n k{Q 2 ) 



2M n M k 



■ J-nk(Q 2 



Q 2 + M 2 + M 2 
For the case of diagonal n — > n transitions this gives 

^ nn{Q ) 



(43) 



F nn (Q z ) 



1 



2 /2M, 



J2 



(44) 



expression similar to that derived in Ref. [18| . 

Thus, we can obtain F nk (Q 2 ) form factors from the 
basic form factors ^F n k(Q 2 )- Note, that these form factors 
also have a generalized VMD representation 



Fnk{Q 2 



oo 



T 

J n 



^ + Q 2 /Ml ' 



(45) 
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with the coefficients T m ^ n k given by the overlap integrals 

9b frr 



m,nk 



M 2 



dzztpm(z) <S> n {z)<S> k (z) 



For the lowest diagonal transition (i.e., for n — k 
we have 



(46) 
0) 



^oo (Q ) = 2 / dzze~ z J(Q,z) 



(47) 



Incorporating the representation lfl4|) for J{Q, z) and us- 
ing a = Q 2 /4k 2 , we obtain 



•foo (Q 2 ) = 



1 



1 



1 



l + Q 2 /M 2 



(48) 



Here, we took into account that the mass of the lowest 
bound state (i.e., p— meson) is M = M p = 2k. 

Notice, that we obtained exact vector meson domi- 
nance for J r oo(Q 2 ): this form factor is completely deter- 
mined by the lowest bound state in the g-channel. The 
higher states do not contribute because the overlap inte- 
gral Frnfio corresponding to the contribution of the m th 
q-channel bound state vanishes for m > 0: 



ro.OO 



dzz 3 e 2,2 L^z 2 ) = 5. 



mO 



(49) 



In the case of F 00 (Q 2 ) form factor, we have 



FaoiQ 2 ) = 



1 



(l + o)(l + a/2) 
2 



1 + Q 2 /M 2 1 + Q 2 /M 2 



(50) 



Thus, the i 7 bo(Q 2 ) form factor is given by contribu- 
tions from the lowest two g-channel bound states. Since 
Foo(Q 2 ) ~ l/<3 4 for large Q 2 , exact VMD is impossible 
for this form factor: other resonances are needed to "con- 
spire" to cancel their leading l/Q 2 terms at large Q 2 . In 
the soft-wall model, this cancellation is provided by just 
the first excited state. 

For small Q 2 , the form factor Fq (Q 2 ) has the following 
expansion: 



F 00 (Q 2 ) = 



2M 2 + AMy W) 



(51) 



The Lorentz structure of the 3-point function in the 
soft-wall model is the same as in the hard-wall model 
considered in Ref. [HI], where it was shown that electric 
Gc, magnetic Gm and quadrupole Gq form factors (for 
definitions, see, e.g., [H, HH) of the n th bound state are 
all expressed through the F nn (Q 2 ) form factor: 



G^{Q 2 )^-F nn {Q 2 ) 

Q 2 



GcHq 2 



6M ; 



g { m\Q 2 ) 



2F nn (Q 2 ) , (52) 



The same relations hold for the soft-wall model. As a 
result, small-Q 2 expansion of the electric form factor of 
the lowest bound state in the soft-wall model is given by 



G 00 (Q 2 ) 



1 



Foo(Q 2 ) 

Q 4 



6M 2 

5W + 2 W + W 



(53) 



and the electric radius for the p-meson in the soft-wall 
model is 



(r 



2\S 
pl 



0.655 fm 



(54) 



This radius is larger than the value (r 2 ) H = 0.53 fm 2 that 
we obtained in Ref. [HI in the case of the hard- wall cutoff. 

The radius of the n th excited state can be found from 
the slope of F nn (Q 2 ). The latter can be calculated using 
Eqs. ([7]), (|3"4"|) . Defining the slope coefficient S n by 



dQ 2 F "" (Q2) 



Q 2 =Q 



M (55) 



and using explicit form of Laguerre polynomials, we find 

c \ ' /Tm+l/ii+1 / 1 \l+rn ( m + [ + jOj 

> ^n+l^n+A i; (n + l)m\l\ 



E 

m,l=0 



m+l+2 



p=l 



P 



(56) 



(Cf are binomial coefficients). A faster algorithm for 
numerical calculations is provided by the formula 



Q \ f-1771 /~1 

rn=0 



n—m 



m-\-n-\-l / ^ ^n — m 
k=0 



2m+fc+2 1 

M) fc - • ( 5? ) 



For n — 0, these expressions give the result So = 3/2 
corresponding to Eq. (|51|) . For higher states, we have 
Si = 23/12, S 2 = 11/5, S 3 m 2.415, S 10 « 3.245, S 20 « 
3.816, S 5Q « 4.633, 5 100 « 5.281, S 150 « 5.667, S 200 « 
5.943 . For n > 2, these values are well approximated by 
a simple empirical formula 



S n «ln(n + l) + ? + - 7 -A 
3 4(n + 



1 



(58) 



Thus, the squared sizes of excited states increase with 
the excitation number n. However, contrary to ex- 
pectations of Ref. [12], the raise is only logarithmic, 
( r n) S ~ Inn rather than linear. Such an outcome is not 
unnatural since Eq. (|56p differs from the identity 



E Cn^M-ir 



m,l=0 



(m + i + iy. 

(n+ l)m\l\ 



= 1 



(59) 



(that follows from the normalization condition ([9])) by 
the sum 



m+l+2 



E 1 



In (m + I + 2) 



(60) 



m-\-l—*oo 
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which has a logarithmic behavior for large m + l, and for 
large n it may be approximated by Inn for the bulk of 
m, I values. However, it would be interesting to derive a 
formal proof. 

It should be noted, that in the hard-wall model , the 
slope of F nn (Q 2 ) at Q 2 — decreases with n. For the 
lowest state, the value S^ 1 = 1.192 was found in Ref. [HI]. 
For higher radial excitations, we have S 1 ^ = 0.877, S% = 
0.833, S? = 0.806, S? = 0.804, S? 00 = 0.803, i.e., {r 2 n ) H 
tends to a constant value as n — > oo. 



VI. COMPARISON WITH HARD- WALL 
MODEL 

Note that in the hard-wall model considered in 
Ref. [l8(, all the g-channel states give nonzero contri- 
butions to JFoo (Q 2 )- In fact, it is strongly dominated by 
two lowest g-channel states. The role of the first exci- 
tation in the hard-wall model is especially important for 
large Q 2 : it gives asymptotically 2.061 M 2 /Q 2 while the 
lowest state contributes only 0.619 M 2 /Q 2 . 

Q 2 F(Q 2 ) 




FIG. 1: Q 2 -multiplied p-meson form factor J-oo{Q 2 ) (dis- 
played in GeV 2 ) as a function of Q 2 (given in GeV 2 ) in hard- 
wall (upper line, red online) and soft-wall (lower line, blue 
online) models. 

It should also be mentioned that in both models 
•7"oo (Q 2 ) has ~ 1/Q 2 behavior for large Q 2 . How- 
ever, the normalization of the asymptotic behavior in 
hard-wall model is much larger than in soft-wall model: 
FoaiQ 2 ) - 2.566 M 2 /Q 2 , while ^ S (Q 2 ) - M 2 /Q 2 . 

As discussed in Refs. [H,[l8j], to calculate the large-Q 2 
behavior of ^^(Q 2 ), one should take the large-Q 2 limit 
of i/ H (Q, z), which is given by zQKi(zQ) = K,(Qz), the 
free-field version of the nonnormalizable mode. Asymp- 
totically, it behaves like e~® z , so only small values of z 
are important in the relevant integral. As a result, 



FmiQ 2 ) 



l<jff(Q)l s 

Q 2 



dxx 2 Kx{ X ) 



2|^(0)| S 

Q 2 



( 61 ) 

i.e., the large-Q 2 behavior of T^ Q {Q 2 ) is determined by 
the value of the <!> wave function at the origin, which is 



given by 



$«(0) 



V2M p 



1.133 M 



(62) 



7o, 1^1(70, 1/ 

The nonnormalizable mode J s (Q,z) of the soft-wall 
model should also convert into IC(Qz) when Q 2 is large. 
To see this directly, we compare the integral representa- 
tion 



lC{Qz) 



dx 



(1-xf 



exp 



{l-x)Q 2 xz 2 



Ax 



1 



(63) 



for IC(Qz) and the representation 



dx 



{l-xf 



exp 



Q 2 



In 



- x 
(64) 



for J s (Q,z) following from Eq. fTi]). For large Q 2 , 
both integrals are dominated by the region where 
1 — x <~ 2zk 2 /Q. Then both (1 — x)/x and ln(l/x) may 
be approximated by (1 — x). Thus, large-Q 2 behavior of 
J s (Q,z) coincides with that of JC(Qz), and Eq. (ffJTj) is 
applicable in soft-wall model as well, with the normal- 
ization of the asymptotically leading term determined by 
the value of $0(2) at the origin, which is 



(65) 



$g(0) = M p /V2 w 0.707 M, 



Hence, it is the difference in the values of <E> wave func- 
tions at the origin that explains the difference in the 
asymptotic normalization of Tqq (Q 2 ) in these two mod- 
els. 

The difference in the values of $(0) leads also to dif- 
ference in the values of coupling constants /„ related to 
*n(0) by 



/„ = *„(0)M n /</ 5 



(66) 



The constant 55 is determined by matching the asymp- 
totic behavior 



^AdS 



P 



(67) 



of the two-point function £ Ads (p 2 ) given by Eq. (I20|) with 
the QCD result for the correlator of the vector currents 
Jfj, = d'y^u having quantum numbers of the p + meson. 
Since 



s QCD (p 2 



12tt ; 



p 2 Hp 2 ) 



we have 



95 



for N c — 3. This gives 



Ml 
2tt 



V2tt 



(309 MeV) 2 



(68) 



(69) 



(70) 
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for the p coupling constant in the soft-wall model, and 



J n 



M 2 



^70,1^1(70,1) 



(392 MeV) 2 



(71) 



in the hard- wall model [28| . The experimental value [22] 



f^P = (401 ±4 MeV) 2 



(72) 



is very close to the hard- wall model result, and in this 
respect the hard-wall model is more successful. It may 
be also noted that, unlike the value (r 2 ) s = 0.655 fm 2 
in Eq. (53]), the hard-wall model result (r 2 ) H = 0.53 fm 2 
for the p-meson charge radius obtained in our paper 
[l8| practically coincides both with the Dyson-Schwinger 
model result of Ref. [23j and lattice gauge calculation re- 
ported in Ref. [24J. 

It is also instructive to consider the modified coupling 
g p = fp/Mp that has the dimension of mass, and deter- 
mines the asymptotical behavior of the form factor. Its 
value in the soft-wall model 



4 



^ w 123 MeV 
2tt 



(73) 



is close to the experimental value of the pion decay con- 
stant fx w 131 MeV. Moreover, the pure p-pole result 
(|48f is close to the experimental data on the pion form 
factor. So, it is tempting to take for the pion the same 
wave functions that were obtained in the p-meson case 
and use Eq. (jig)) as a model for the pion form factor. 
This was done in the paper |25j (that appeared after we 
submitted the original version [26j of the present paper 
to the arxive). Taking k = 375 MeV (which is slightly 
smaller than m p /2), the authors obtained good agree- 
ment of the 1/(1 + Q 2 /4k 2 ) curve with the pion form 
factor data (though the value of f 2 is then about 30% 
below the experimental one). However, within the model 
of Refs. [f| [a [Hj], which we follow here, the analysis of 
the axial- vector current channel requires the inclusion of 
chiral symmetry breaking effects absent in the vector cur- 
rent channel. As a result, wave function equations for the 
pion are completely different from those for the p-meson. 
We discuss the pion form factor in a separate publication 

M. 



VII. SUMMARY 

In the present paper, we studied wave functions and 
form factors of vector mesons within the framework of the 
soft-wall model [HJ which produces a more realistic spec- 
trum for higher excited mesons [19( than the hard-wall 
model of Refs. [y, [7J, |8( . Our analysis uses the approach 
similar to that we developed in Ref. [18] in application 
to the hard-wall model. 

An essential element of our study of the soft-wall model 
is the integral representation, which we found for the 
bulk-to-boundary propagator V(p, z). It allows to write 



V(p, z) as an explicit expansion over bound state poles. 
In this sense, it plays the same role as the Kneser- 
Sommerfeld expansion that we used in our study [18] of 
the hard-wall model. 

The pole expansion of V(p, z) involves ll ip wave func- 
tions" that describe ^-dependence of a particular pole 
contribution and are eigenfunctions of the 5D equation 
of motion. However, since ij) n {z) wave functions are not 
direct analogues of the usual quantum-mechanical wave 
functions, we introduced wave functions" resembling 
wave functions of oscillator states in quantum mechan- 
ics. In particular, the values of these functions at the 
origin give the couplings g^fn/Mn of the bound states, 
and these functions exponentially decrease with z 2 . 

Analyzing the three-point function, we obtained ex- 
pressions for transition form factors both in terms of the 
ijj wave functions and the "more physical" $ wave func- 
tions. We demonstrated that, just like in the hard-wall 
model, the form factors can be written in the form of 
generalized vector meson dominance representation, i.e., 
as a sum over all bound states in the channel of elec- 
tromagnetic current (this result confirms the claim [l6| 
that generalized VMD is a common feature of AdS/QCD 
models). 

We derived an explicit expression for p-meson form fac- 
tors, and analyzed their behavior in the regions of small 
and large Q 2 . In particular, we calculated the p-meson 
electric radius in the soft-wall model, and found that it 
is larger than in the hard-wall model (the latter agrees 
with calculations in Dyson-Schwinger model [231 ] and lat- 
tice QCD [13]). Our calculation also demonstrated that 
the squared radii of higher excited states increase with 
n, the number of the radially excited level. However, 
contrary to expectations of Ref. [e|] , the increase is only 
logarithmic rather than linear. Another result is that, 
in the soft- wall model, the p-meson form factor T P (Q 2 ) 
(corresponding to the overlap of the <& wave functions) 
exhibits an exact VMD pattern, i.e., it is given by a single 
monopole term due to the lowest bound state in the Q 2 - 
channel. In the case of the p-meson form factor F p (Q 2 ) 
(that is given by the overlap of the ^-wave functions), 
we found a two-resonance dominance pattern, when just 
two lowest bound states in the Q 2 -channel contribute. 

Analyzing the large-Q 2 behavior of the J r p (Q 2 ) form 
factor (given by exact p-pole VMD), we established that 
its asymptotic normalization in the soft-wall model is 
much lower (by factor 2.566) than that of the hard-wall 
model. This difference is explained by essentially lower 
value of the soft- wall model $ wave function at the origin. 

Finally, we calculated the p-meson coupling constant 
f p both in the soft-wall and hard-wall models, and found 
that the experimental value is closer to the hard-wall 
model result. 
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The hard-wall model result F, 



1/2 



329 MeV presented 



in Ref. [6j corresponds to the (u7 M it — d r y ti d)/2 current 
which differs by from the current (u7 M u — dj^d)/^ 
that has the same normalization as d'y^u. 



